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<—| \ Abstract 

A feedback stabilization scheme to stabilize a classical reacting Hamiltonian sys- 
tem is proposed. It is based on transforming a saddle-type equilibrium to an asymp- 
totically stable one, and is given in a simple and algorithmic way. The question of 
destabilization of a stable system to make a reacting system is also addressed. The 
theory is illustrated with the examples of a model Hamiltonian of the form kinetic 
plus potential, and the hydrogen atom in crossed and magnetic fields. 

1 Introduction 

Feedback stabilization of nonlinear systems is a well-established topic in control the- 
ory pp. As a special case the Hamiltonian stabilization appears as a nice technique 
due to the rich geometric structure behind it: Poisson structures. A detailed formula- 
tion of Hamiltonian stabilization can be found in [2] with an emphasis to mechanical 
systems. On the other hand, a more general treatment which extends the stabi- 
lization method to Poisson manifolds is introduced in [3], and this generality led an 
application to systems with symmetry. For a related work on stabilization of me- 
chanical systems with symmetry we refer to [H [5] . A more recent work [6] outlines 
controlling dissipation- induced instabilities. 

Reaction-type dynamics has had a renewal of understanding after the develop- 
ment of its phase space geometric picture [7]. It is based on identifying geometric 
structures which govern reaction dynamics around a saddle-type equilibrium. Since 
the introduction of these structures there have been a big log of work done on the 
dynamics of these systems which are reviewed in [8] . These systems not only include 
chemically reacting systems but also systems in celestial mechanics, atomic physics, 
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diffusion dynamics in materials, which have reaction-type dynamics [8]. On the other 
hand any attempt to stabilize these systems is apparently lacking. 

Our aim in this paper is to address the problem of stabilization of reacting systems 
by means of Hamiltonian stabilization tools. By the special character of reacting sys- 
tems, the problem turns out to be making a saddle-type equilibrium asymptotically 
stable by adding some suitable feedback. The feedback are algorithmically derived 
in terms of linearization of the Hamiltonian vector field. We also touch upon the 
problem of the other way around, namely, one can adopt what is done for saddle 
equilibria to make a center-type equilibria saddle-type. Two examples elucidate the 
theoretic part of the paper: a simple reacting system with a model potential, and 
the hydrogen atom in crossed and magnetic fields. 

We structured the paper as follows. Sec. [2] briefly recalls Hamiltonian stabilization 
and Sect. Ogives a short introduction to geometric theory of reactions. We give a 
detailed explanation of stabilization of reacting systems in Sec. 01 and a brief look at 
destabilization of stable systems in Sec. [5j Examples are given in Sec. [6] which are 
succeeded by conclusions. 

2 Hamiltonian stabilization 

We begin with a brief overview of some basics of the theory of Hamiltonian stabi- 
lization on Poisson manifolds here. A more detailed information can be found in 

mm- 

Let P be Poisson manifold and let H : P — > M. be a Hamiltonian function with 
the corresponding Hamiltonian vector field Xh- Then the equations of motion read 
z = Xh{z). If zo is an equilibrium, i.e. DH{zq) = 0, the Hessian D 2 H(zq) is 
intrinsically defined. Now adding some inputs Fi : P — >■ K, i = l,...,m, such that 
Fi(z ) = 0by 

z = X H (z) - X Fl {z)ui X Fm (z)u m (1) 

gives a system of which zq is an equilibrium as well. Here the feedback are assumed 
to be Ui = di X Fi , i = 1, . . . ,m, for scalars di > 0. Note that, if Fi are in involution, 
i.e. {Fi, Fj} = for all i,j = l,...,m, then Ui = di Fi. 

Associated to the closed- loop system given above the following space is defined: 

C = span^, {H, Fi}, {H, {H, Fi}}, ...,}, (2) 

i = l,...,m. Here the coefficients of the linear combinations are real numbers. 
Accordingly, one defines the associated co-distribution 

dC = span{dg(z)\g £ C}. (3) 

Then we recall the following key result [2] for our application. 

Theorem 1. Let D 2 H(zq) be positive definite and dim(dC) = 2n around zq. Then 
the feedback Ui = Ci X Fi makes zq an asymptotically stable equilibrium. 

The idea behind the proof is to use LaSalle's Principle where the Lyapunov func- 
tion is assumed to be the Hamiltonian, and the dimensionality condition makes sure 
that the only trajectories that lie in a certain neighborhood are the equilibria [9]. 
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Remark 1. A system for which D 2 H(zq) is not positive can also be stabilized under 
suitable conditions Jj| But we want to make D 2 H(zq) positive definite practically 
by using the linearization below. 

Remark 2. In general, it is not easy to show whether dC(z) is constant dimensional 
or 2n- dimensional. But if Fi are independent, m = n, and in the form kinetic plus 
potential, then dim(dC) = 2n is guarantied JIJ/. 

3 Geometry of reacting systems 

In this section, we outline the geometric theory of reaction dynamics briefly as in the 
form given in [10]. A detailed explanation can be found in 0[8], for instance. 

3.1 The linear case 

Consider the simplest reaction-type Hamiltonian , i.e. the quadratic Hamiltonian 
given by 

%?)4^-^+Et(^^- ( 4 ) 

k=2 

where \,0Jt > 0. Then DH(0) = and the matrix associated with the linear vector 
field has real eigenvalues ±A and complex conjugate imaginary eigenvalues ±iwfc, 
k = 2, ... ,n. Integrability of the system can be seen by the constants of motion 

2a = Pi ~ q\ , Zk=pl + ql, k = 2,...,n. (5) 

Consider a fixed energy E > 0, where is the energy of the saddle. Setting q\ = 
on the energy surface gives the (2n — 2)-dimensional sphere 

Sit 2 = {(q,p) € R 2n : H 2 (q,p) =E, qi = 0}. (6) 

The dividing surface S^~ 2 divides the energy surface into the two components which 
have qi < (the 'reactants') and q\ > (the 'products'), respectively, and as q\ = 
dH2/dpi = Xpi 7^ for p\ ^ the dividing surface is everywhere transverse to the 
Hamiltonian flow except for the submanifold where q± = p\ = 0. For q\ = p\ = 0, 
one obtains Y^Jk=2 ^t(Pk + 9fc) = The submanifold thus is a (2n — 3)-dimensional 
sphere which we denote by 

Snhi~m = {(q,p)eM. 2n : H 2 (q,p) = E,q 1 =p 1 =0}. (7) 

This is a so called normally hyperbolic invariant manifold [llj (NHIM for short), i.e. 
^nhim ^ s invariant (since q\ = pi = implies qi = p\ = 0) and the contraction and 
expansion rates for motions on are dominated by those components related 

to directions transverse to fff ^j^ . The NHIM (J7|) can be considered to form the 
equator of the dividing surface © in the sense that it divides it into two hemispheres 
which topologically are (2n — 2)-dimensional balls. All forward reactive trajectories 
(i.e. trajectories moving from reactants to products) cross one of these hemispheres, 
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and all backward reactive trajectories (i.e. trajectories moving from products to 
reactants) cross the other of these hemispheres. Note that a trajectory is reactive 
only if it has X\ > (i.e. if it has sufficient energy in the first degree of freedom). 
Trajectories with X\ < are nonreactive, i.e. they stay on the side of reactants or 
on the side of products. See Fig. [1] for the phase portrait of the system. A forward 
reactive trajectory is depicted by the dashed curves. 



3.2 The general (nonlinear) case 

Consider a Hamiltonian H : R 2n — > R with an equilibrium at the origin (q,p) = (0, 0) 
for some canonical coordinates (q,p). Assume that H has a saddle-center-. . .-center 
stability type equilibrium, that is, the matrix associated with the linearization at 
(0,0) of the Hamiltonian vector field has eigenvalues =fA, A > 0, and ^ia;^, w& > 0, 
k = 1, . . . ,n — 1. Further assume that the submatrix corresponding to the imaginary 
eigenvalues is semisimple. In the neighborhood of the saddle the dynamics is thus 
similar to that of the linear vector field described in Sec. 13.11 In fact if follows from 
general principles that all the phase structures discussed in Sec. 13.11 persist in the 
neighborhood of the saddle (which in particular implies that one has to restrict to 
energies close to the energy of the saddle). Moreover, these phase space structures 
can be constructed in an algorithmic fashion using a Poincare-Birkhoff normal form 
[8]. Assuming that the eigenvalues LUk, k = 2, . . . , n, are independent over the field 
of rational numbers (i.e. in the absence of resonances), the Poincare-Birkhoff normal 
form yields a symplectic transformation to new (normal form) coordinates such that 
the transformed Hamiltonian function truncated at order no of its Taylor expansion 
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assumes the form 

Hnf(Xi,I2, . . . ,Z n ) , (8) 

where X\ and Ik, k = 2, ... ,n, are constants of motions which (when expressed in 
terms of the normal form coordinates) have the same form as in ([5]), and H^f is a 
polynomial of order no/2 in X\ and I k , k = 2, . . . , n (note that only even orders no 
of a normal form make sense) . 

In terms of the normal form coordinates the phase space structures can be defined 
in a manner which is virtually identical to the linear case by replacing H2(q,p) by 
H nf (Ii, I2, ■ ■ ■ ,X n ) in the definitions in Sec. 13,11 Using then the inverse of the normal 
form transformation allows one to construct the phase space structures in the original 
('physical') coordinates. As it is seen the Poincare-Birkhoff normal form is the main 
tool in defining phase space structures, but we only review the first order linearization 
in the next section, as it serves enough for our purpose of stabilization. 



4 Stabilization of saddle-type equilibria 

In this section we consider a nonlinear Hamiltonian system of reaction-type, and we 
give a result which gives an algorithmic way of making the system asymptotically 
stable around the given equilibrium. This is done as follows. 

Let P be a Poisson manifold, (z±, . . . , z 2n ) = (xi, ■ ■ ■ , x n , Pi, ... , P n ) be a canon- 
ical coordinate system around zq £ P which is set to be the origin (0, ... ,0), and 
H : P — > R be an analytic Hamiltonian function with the corresponding Hamiltonian 
vector field Xn having zq as an equilibrium point of type saddle-center-. . . -center. So 
we assume that the linearization matrix of Xn, or in other words the matrix J D 2 H 
evaluated at zq, has eigenvalues =fA, ^iuJk, for reals A, uj^ > 0, k = 2, . . . , n. 

One can put the quadratic part of H into the form Q by a symplectic change of 
coordinates [5J. To do this label the eigenvalues by 

ei = X = -e n+ i, e k = iuj k = -e k+n , k = 2,...,n, (9) 

and corresponding eigenvalues by v\, . . . , V2 n - Consider the following symplectic ma- 
trix 

M = (ci v 1 ,c 2 Rev 2 , ■ ■ ■ ,c n Rev n , c\ vi +n , c 2 Imu 2 , . . . , c n \mv n ) (10) 

where 

ci = {vi, Jv 1+n y 1/2 , c k = (Rev k , J\mv k )~ l/2 , k = 2,...,n. (11) 
Then the coordinate transformation 



X 




Q 


p 







giv6S ci new Ccinoniccil coordinate system . . . , Pi,. . . p n ), and in these coordi- 
nates the quadratic part of the Hamiltonian takes the form 

n 

H 2 {q,p) = \q 1 p 1 + Y,'^(fk + ?k)- (13) 

k=2 
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Let a rotation N be given by 

1 



1 



then if we introduce 
the coordinate transformation 



qk = Qk, Pk=Pk, k = 2,...n, 
S = N ■ M~ l 



S- 



X 




Q 


p 







(14) 

(15) 
(16) 



gives also a set of canonical coordinates (q\ . . . ,q n ,p% . . . ,p n )- Then in (q,p) coordi- 
nates H2 assumes the form (0]). 
Consider the controls 



-cFi — diFi, i = l,...,n, (17) 

, n, where the functions 

(18) 



with any constants such that c > A and di > 0, i = 1, 
Fi : P — > R are given by 

n 

Fi = $ij z j 

for the matrix S introduced in Eq. [15j Then we prove 
Theorem 2. With the notion above, the system 

z = X H (z) - X Fl (z)ui X Fn (z) u r 

is asymptotically stable around zq. 



(19) 



Proof. We want to show that the conditions of Theorem [T] are satisfied. First, it 
is easily seen that Fi(zo) = Y^j=i &ij z j( z o) = for i = 1, . . . , n. Then we add the 
control v = —cF\ and consider the new system 



z = X H - X Fl (z) v. 



(20) 



One can check that the system (|20p is also Hamiltonian with the modified Hamilto- 
nian i^ moc j = H + \cF^ and the Hessian matrix D 2 H(zq) is positive definite. In 
fact, 

z = X H - X Fl (z)v = X H + c X Fl (z) F 1 =X H + cXx_ F? (z) = X H+1 _ c F? (z) , (21) 
and in normal form coordinates F{ = i = 1, . . . , n, so in these coordinates 

^H mod (q,p) = d2h ^p) +cE 11} (22) 

where En is the matrix with zero entries except the first entry equal to 1. This 
shows that 

D2H mod( z o) = diag(-A + c,\,ui,Ui, . . .,uj n ,uj n ). (23) 
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which is, clearly, positive definite since c > A. 

So far, we have ensured the positive definiteness condition in Theorem[TJ Next we 
add the controls Vi = —di F(, i = 1, . . . , n, and it remains to show that the functions 
Fi, i = 1, . . . , n satisfy the dimensionality assumption dim(dC) = 2n for 

dC(z) = span{dg(z)\g G C}, (24) 

where 

C = span{F i ,{H mod ,F l },{H,{F t }},...}, i = l,...n, (25) 
and Fi are in involution. To see this it is observed that 

C = span{g i ,{ J ff mod oS' _1 ,ft},{i?,{ft,}},...}, i = l,...n. (26) 

In fact, by Fi = q^o S we have 

i H mod> F i} = { H mod> H oS } = { H mod ° ?i) ( 27 ) 

since 5 is a Poisson map. But -ff moc ; ° <5 _1 is real analytic so we can write it as a 
Taylor series around zq where the quadratic part is given by (H]). Then it can be seen 
that fi = {H mQ( j L o S^ 1 , qi}, i = 1, . . . , n, are independent, because one has 

h = {#mod °S-\q 1 } = \ Pl + h.o.t., 

fk = {H mod o S^^k} =u k p k + h.o.t., k = l,...,n. 

Furthermore, {qi, . . . , q n , f±, . . . , f n } forms a set of 2n independent functions. So, 
dC(z) is 2n-dimensional, in particular 

dC(z ) = span{^i(z ), . . . ,dq n (z ),dp 1 (z ), . . . ,dp n (z )}. (29) 

As the final step we need to show that Fi are in involution. This can also be seen 
easily by 

{F^Fj} = { qi o S, Qj o S} = {q h qj } = 0. (30) 

□ 

The feedback added system is no more conservative because of the dissipative 
inputs Vi. As the system is asymptotically stable, trajectories projected into the 
phase planes look like the ones in Fig. [2j 

5 Destabilization of a stable system to make a 
reacting one 

A similar procedure as in Sec. H] can also be applied to a stable system with purely 
complex eigenvalues in order to obtain an unstable system with a saddle. 

We cansider again a Poisson manifold P, a canonical coordinate system denoted 
by (z±, . . . , Z2n) = ■ ■ ■ , x n , Pi, ... , P n ) around zq G P which is set to be the 
origin (0, . . . , 0), and be an analytic Hamiltonian H : P — )• 1R with the corresponding 



Figure 2: Projetion of trajectories near the equilibrium into the phase planes. 



Hamiltonian vector field Xh having zq as an equilibrium point of type center-. . . - 
center. So we assume that the linearization matrix of Xjj, or in other words the 
matrix J D 2 H evaluated at zq, has eigenvalues ^iwj, for reals > 0, i = 1, . . . , n. 
The quadratic part of the Hamiltonian can be put in the form H into the form 



i=l 



(31) 



by a symplectic change of coordinates [8] as recalled in the following. Label the 
eigenvalues by 

ei = iwj = -e i+n , i = l,...,n, (32) 

and corresponding eigenvalues by v\, . . . , V2 n - Consider the following symplectic ma- 
trix 

M = (ci Kev\, . . . , c n Rev n , c\ Imt>i, . . . , c n lm.v n ) (33) 



where 

Ci = (Revi, Jlmvi}~ 1/2 , i = l,..., 
Then the coordinate transformation 

M ■ 



n. (34) 



(35) 



gives a new canonical coordinate system (qi, . . . , q n ,pi, ■ ■ ■ p n ), and in these coordi- 
nates the quadratic part of the Hamiltonian takes the form (|3ip . 
Consider the control 

u = -cF 1 (36) 



X 




Q 


p 




. P . 
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with any constants such that c > uj\ the function F\ : P — > R is given by 

n 

F 1 =Y,MijZ j (37) 

for the matrix M introduced in Eq. [35j Then we prove 
Theorem 3. With the notion above, the system 

z = X H {z)-X Fl {z)u 1 (38) 
is of type saddle — center — ■ ■ ■ — center. 

Proof. Similar to the proof of Theorem [2j □ 

6 Examples 

We illustrate the procedure of making a saddle- type equilibrium asymptotically stable 
with two examples. 

6.1 A model example 

The following system represents a typical pattern of isomerization reactions [12J. 
Consider a system with potential function 

V(xi,x 2 ) = -$x\ {x\ - l) 2 + jkx\, a>b, (39) 
a z o z 

and Hamiltonian 

H{x x , x 2 , Pi, P 2 ) = ^ Pi + ± Pi + V( Xl ,x 2 ). (40) 

Clearly, H has three equilibria which are critical points of V. These are two centers; 
(xi,X2) = (0,0) and (x\,X2) = (1,0), and one saddle; (xi,x 2 ) = (1/2,0). As we are 
interested in the saddle, we translate the coordinates by (xi,x 2 ) i-> (x\ + 1/2, x 2 ) 
to move the saddle to the origin. We use the same notation for the translated 
coordinates and the potential, then we have 

V{xi,x 2 ) = \ x\(xi - l) 2 + ^x\, a>b (41) 
a z b z 

which will be used henceforth. The contours of the potential surface are depicted in 

Fig. m 

The Hessian matrix J D 2 H at point z$ = (0,0,0,0) is computed to be 
JD 2 H(z ) = 







i 

771 







2 

"7? 



1 

1 







(42) 
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Figure 3: Contours of potential functions; (a) the original model potential (HT]1 . (b) the 
modified potential 



which has eigenvalues e\ = ^ = — e 3 , e 2 = -^p i = — e4 with corresponding eigen- 
vectors v x = {a, 0, 1, 0}, v 3 = {-a, 0, 1,0}, v 2 = {0, i, 0, 1}, v 4 = {0, i, 0, 1}. 

So, we have cj = {v x ■ Jv 3 y 1/2 = (2a)- 1 / 2 , c 2 = (Re(u 2 ) • Jlm^))" 1 / 2 = 2 1 /H- 1 /2. 
Then the matrix M reads 



M 




1 

y/2a 










21/4 

>/6 



v» 





1 

/2a 




V6 
2^74 






(43) 



Finally, the matrix S which is the multiplication of the rotation matrix 



R 



1 



1 





and the matrix M becomes 



S 








2V 4 

' 



1 


1 

















Vb 




(44) 







Hence the functions F±, F 2 are obtained to be F\ 
the modified Hamiltonian -ff moc | has the form H mo( ^ 



XI 



F 2 = 



Vb 
2 1 / 4 



1- 



(45) 



P 2 . Observe that, 
and the modified 
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potential is 



V, 



mod 



V + 



2 a 2 



T 2 

x 1 



(46) 



of which contours are depicted in Fig. [3](b). After the addition of associated controls, 
we have the system with equations of motion 



where we choose d\ = d 2 = 1. 



X 2 

A 

P-2 



P 



P2 + 



V2 



X2 



4x? Pi 
H 

a 



a* 
2x 2 

' b 2 ' 



(47) 



6.2 Hydrogen atom in crossed and magnetic fields 

The following example is a Hamiltonian system which is not of the form kinetic 
plus potential. We do not give the original form but a form obtained after some 
manipulations [7J. 

The Hamiltonian can be put in the form 



h= l - (A 2 + p| + p 3 2 ) - 1 + J (xiP 2 -x 2 p 1 ) + J (xf + 4) 



exi, 



(48) 



where R 



(xl 2 + xl + x|) 1/2 . 



We will consider the experimentally interesting 



value e = 0.58 henceforth. The Stark saddle point in atomic physics corresponds to 
the point (xi, x 2 , x%, Pi, P 2 , P3) = (e 1 / 2 , 0, 0, 0, — \ e -1 / 2 ). So after a coordinate shift 
(x\ , x 2 , x 3 , Pi , P 2 , P 3 ) (->■ (xi - e 1/2 , x 2 , x 3 , Pi , P 2 , P 3 + \ er 1 / 2 ) , by retaining the same 
notation for the translated coordinates and the new Hamiltonian we have 



// - (A 2 + P| + P| ) " I + I (*i P2 - x 2 Pi) + ^ (x 2 + x 2 ) 



=1/2 



(49) 



1 /2 

where P = [(x\ + e 1//2 ) 2 + x 2 + x 2 ] . Then the matrix J D 2 H(0) is obtained to be 



JD 2 H(0) 





0.5 


0.63343 





-0.5 




-0.691715 










-0.441715 



1 




0.5 






1 
1 

-0.5 







(50) 



which has eigenvalues 

e x = 0.63645 = -e 4 , e 2 = 0.981506 i = -e 5 , e 3 = 0.664616 i = -e 6 (51) 
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with corresponding eigenvectors 



So, we have 



VlA = {±0.63645, 0.478361, 0, 0.644249, ±0.0137719, 0}, 
v 2 , 5 = {±0.981506 i, 1.84678, 0,-0.0399619, ±1.32188 i,0}, 
^3,6 = {0, 0, ±1.50463 i, 0,0,1}. 

Cl = (v! ■ Jv A )- l/2 = 1.09551, 

c 2 = (Re(v 2 ) ■ JIm(t> 2 ))~ 1/2 = 0.634944 

c 3 = (Re(u 3 ) • JIm(v 3 ))~ 1/2 = 0.81524. 



(52) 



(53) 



M 



Then the matrix M reads 

0.697235 -0.697235 0.623201 

0.524048 1.1726 0.524048 

-1.22663 

0.705779 -0.0253736 0.705779 

-0.0150872 0.0150872 0.839317 

0.81524 

(54) 

Finally, the matrix S which is the multiplication of the rotation matrix R and the 
matrix M becomes 



0.998122 




0.0253736 






0.839317 


0.0213366 











-0.81524 





-0.623201 


0.986039 





-0.741116 




1.1726 







1.22663 






Hence the functions F± , F2 , F3 are derived as 

F x = -0.741116 P 2 ±0.998122xi, 
F 2 = -0.623201 PI ± 0.839317x2, 
F 3 = 1.22663 P 3 . 

This way, the new system is made asymptotically stable around the origin. 



(55) 



(56) 



7 Conclusions and future work 

An algorithmic stabilization of reacting systems is outlined. It relays on the lin- 
earization of the Hamiltonian vector field around the equilibrium. The examples 
reflect the novelty of the technique given in the paper. Next step is to do a study 
for Hamiltonian systems with symmetry where the equilibria are replaced by relative 
equilibria. This can be done by using canonical coordinates on the reduced space 
instead of the reduced energy momentum method as in [3]. A derivation method of 
canonical coordinates on a reduced space for A^-body reduction is outlined in [TO] 
and for cotangent bundle reduction is given in [13j . 
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